Approximate scaling relation for the anharmonic electron-phonon problem 
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An approximate scaling relation is found for the transition 
temperature to a charge-density-wave instability in the anhar- 
monic electron-phonon problem, which maps a wide range of 
interaction strengths, anharmonicities, and phonon frequen- 
cies onto a common functional form. The relation employs the 
wave- function renormalization parameter and is valid even for 
systems that are not Fermi liquids. 



Introduction The interaction of electrons with anhar- 
monic lattice vibrations is a long-standing problem that 
is not yet fully understood. What is surprising about this 
problem is that nearly all real materials are anharmonic 
(as can be seen by the fact that they expand or con- 
tract upon heating), but quasiharmonic models (which 
replace the anharmonic phonons by harmonic phonons 
with temperature-dependent phonon frequencies) work 
remarkably well at describing properties of most materi- 
als 01 . Superconductivity is described most accurately, 
where the theory of electrons interacting with harmonic 
phonons, introduced by Migdal |g] and Eliashberg [g[, 
can routinely reproduce experimental tunneling conduc- 
tances to better than one part in a thousand. The ex- 
planation for this result is actually quite simple — the 
thermal effects that arise due to a nonuniform spac- 
ing of the anharmonic energy levels, are unimportant 
when the temperature is much less than the effective 
phonon frequency (defined by the difference in energy 
between the ground and the first-excited state of the 
anharmonic phonon) 0]. Furthermore, quantum Monte 
Carlo (QMC) studies B, have shown that anharmonic- 
ity does not appear to produce any exotic behavior, such 
as enhancements of transition temperatures, or novel su- 
perconducting behavior Q . Instead, the results indicate 
that the main effect of anharmonicity is to generically 
break particle-hole symmetry. The discovery we present 
here is that the anharmonic systems can be mapped onto 
harmonic ones, with results from widely different param- 
eter regimes collapsing onto the same scaling curve. We 
believe that this result sheds light onto the question of 



why harmonic models work so well for describing prop- 
erties of real materials. 

Our strategy is to solve anharmonic models in the limit 
of large spatial dimension |H] where the lattice many- 
body problem can be mapped onto a self-consistently 
embedded impurity problem that is solved via a QMC 
simulation |^ for quantum phonons, or via an iterative 
transcendental equation for classical phonons |^,|l^ • 

Model The simplest electron-phonon model that in- 
cludes anharmonic effects is the anharmonic Holstein 
model |ll[| in which the conduction electrons interact 
with local phonon modes: 



H = -Y^ Ujcl^Cja + ^igxi - ti){ni^ + n,x) 






(1) 



Here, c]^ (cia) creates (destroys) an electron at site i with 
spin (7, riia- = Ci^Cia is the electron number operator, fl 
is the chemical potential, and Xi (pi) is the phonon coor- 
dinate (momentum) at site i. The hopping of electrons 
is restricted to nearest-neighbor lattice sites i and j on 
a d-dimensional hypcrcubic lattice, tij is isotropic with 
magnitude t =: t* /2\/d where i* = 1 to define the energy 
scale. The bare density of states then becomes a Gaus- 
sian exp(— e^)/Y^. The local phonon has a mass M and 
a spring constant R, associated with it. The anharmonic 
contribution to the phonon potential energy is chosen to 
be a quartic in the phonon coordinate with a strength 
aan- The deformation potential (electron-phonon inter- 
action strength) is parameterized by an energy per unit 
length g. 

A shift of the phonon coordinate is useful for calcula- 
tions, and for illustrating the particle-hole symmetry of 
the model. We shift Xi =: Xi+x' , with g+nx' +3aanx'^ — 
0, to transform the Hamiltonian into 

H = -Y^ Ujcl^Cja + ^{gxr - A*)("«T + "a - 1) 






2M 



(2) 



^ 1999 by the authors. Reproduction of this article by any 
means is permitted for non-commercial purposes. 



with K :~ R,+ I2aanx'^i Pan '■— 'ictanx' , and ^ :~ fl — gx'. 
It is the presence of the cubic term when aan 7^ that 
removes particle-hole symmetry from the problem when 
/J, = 0, since the particle-hole transformation is Xi — > —Xi 



and riia ^- 1 — n^cr- The system does continue to pos- 
sess nesting at half-filling, though, which implies that it 
will have a nonzero transition temperature to a charge- 
density-wave (CDW) at half-filling for any nonzero inter- 
action strength g. 

We examine two different cases in this contribution: 
quantum phonons with a large enough phonon frequency 
that vertex corrections are important (f2 = \J kJM = 0.5, 
corresponding to approximately one eighth the effective 
bandwidth) and classical phonons with zero phonon fre- 
quency [M —f oo and f2 = 0). The former problem is 
solved using QMC techniques that have been described 
elsewhere ||5|,|8|. Transition temperatures are determined 
by calculating the relevant susceptibility (in this case to 
a chessboard-phase CDW) , and determining the temper- 
ature where it diverges (when Trotter error is important, 
we extrapolate Tc using At = 0.2 and 0.4, otherwise, 
we use the larger At). In the latter case, the problem 
is solved within the ordered phase, employing a gener- 
alization of the Brandt-Mielsch formalism to the static 
Holstein model |§,|I^ , and determining Tc as the highest 
temperature that sustains long-range order. 

Results We begin in Fig. f(a) by showing the tran- 
sition temperature to the commensurate CDW at half 
filling for the harmonic case, and three different phonon 
frequencies, f2 = 0, f2 = 0.5, and fi = cx) which is 
identical to the attractive Hubbard model [ppj]. The 
symbols are the exact, or QMC results (the dotted 
lines are guides to the eye), while the other lines are 
second-order weak-coupling conserving approximations 
[ fi3[ (dashed) and second-order strong-coupling calcula- 
tions iQ (solid). Notice how the maximal Tc is essen- 
tially independent of phonon frequency, and that the 
phase diagrams are not too sensitive to phonon frequency 
when il is smaller than the bandwidth. However, note 
that the classical phonon case {fl — 0), and the quan- 
tum phonon case are quite different, since T is always 
above fl for the classical phonons, but T < il for the 
quantum phonon cases shown here. The similarity in Tc 
is surprising, because the CDW vertex is strongly tem- 
perature dependent for the quantum phonons, with it's 
magnitude changing between —Ig^ j k. and —g^ / k as the 
temperature and frequency are varied. 

The anharmonicity is turned on in Fig. f (b) for similar 
effective electron-electron coupling strengths and il = 
and f2 = 0.5. The infinite-phonon frequency limit is 
unaffected by the anharmonicity, since the phonons re- 
spond instantaneously to the electrons, and remain at 
the origin. We plot the strength of the anharmonicity as 
ttonS*^, where x* is the equilibrium coordinate of the 
atomic system with one electron per site, found from 
Eq. (nl) with tij = 0. In the case of weak anharmonicity, 
this equilibrium coordinate lies near —g/n and the anhar- 
monic energy grows linearly with aam but as the anhar- 
monicity becomes large, then x* lies near —{g/ctanY, 



and the anharmonic energy decreases as aan ■ This is 
why the curves in Fig. (f (b) are multivalued for some val- 
ues of the interaction strength. Notice that the results 
for different phonon frequencies are similar in qualitative 
behavior, but that the quantitative results can differ by 
large amounts as the anharmonicity increases. 

We find that the Tc satisfies an approximate scaling 
relation when we plot it instead as a function of 1/Z(0), 
the wave-function renormalization parameter, in Fig. 2. 
Nearly two-hundred data points collapse onto the same 
scaling curve. But the scaling is only approximate, be- 
cause the infinite-phonon frequency limit, lies well off the 
scaling curve for low-to-moderate phonon frequencies, so 
this result must break down as fi increases (note, how- 
ever, that the value fi/i* — 0.5 is larger than the phonon 
frequency in nearly all real materials). The other lines 
in Fig. 2 correspond to different weak-coupling approxi- 
mations plotted now as a function of Z{Q) rather than of 
U . The wave function renormalization parameter Z{Q) is 
extracted from the calculations by a linear extrapolation 
along the Matsubara frequency axis. We compute 

3 ImS(ii7rr) l ImS](3z7rr) 
^^^>-^-2 ttT +2 3^r ' ^^> 

and evaluate it at a temperature T just above Tc- This 
imaginary-axis extrapolation procedure is robust in pro- 
ducing a scaling result even if the system is not a Fermi- 
liquid, where 1/Z(0) would measure the quasiparticle 
weight. For example, in the classical-phonon case, the 
lowest-order contribution to the self energy is S(z) — 
—TUG(z). The curvature oiImS{uj) at the chemical po- 
tential is positive. Hence the system is never a Fermi 
liquid, except at T = 0, where it becomes noninter- 
acting. Nevertheless, the extrapolation procedure given 
in Eq. (g) still falls on the scaling curve, and for the 
weakly-coupled classical-phonon case, it produces Z{0) = 
1 +4:\U\p{0)/3, rather than the expected result of Z{0) = 
1 -I- |C/|p(0) from Migdal-Eliashberg theory. We find that 
a phenomenological functional form that fits the data rel- 
atively well is Tc = O.f 82Z(0)-°i2 exp[-0.75/{Z(0) - 1}] 
(the solid line in Fig. 2). 

As a consequence of this scaling, if we use the quasihar- 
monic approximation to describe an anharmonic system 
with the parameters adjusted to obtain correspondence 
in Z(0), then we also accurately reproduce two-particle 
properties like Tc- Since this result is robust against in- 
troducing anharmonicity, we believe this is the reason 
why the quasiharmonic approximation works so well in 
real materials. 

We plot the interacting density of states and the self 
energy in Figs. 3-5 for weak-coupling harmonic and an- 
harmonic cases, and for a strong-coupling harmonic case. 
We range from a high temperature down to just above 
Tc- Notice how the weak-coupling cases are not Fermi 
liquids, and have quite different DOS and self energies, 



but both map onto the same Z{0) and Tc [our defini- 
tion of Z(0) from the imaginary axis is not equal to the 
derivative of the self energy on the real axis, which has 
the opposite sign for a Fermi liquid here]. Furthermore, 
the paramagnetic phase of this model becomes noninter- 
acting as T — !■ since S vanishes in the weak-coupling 
limit, which also explains why the anharmonic DOS be- 
comes more symmetric as T ^ 0. Scaling even holds (in 
a more approximate way) when the system is a bipola- 
ronic insulator as shown in Fig. 5 — the self-energy has a 
large imaginary part near the chemical potential, which 
grows as T is lowered, and a pseudogap develops in the 
interacting density of states. This is then followed by a 
gap developing at lower T and the imaginary part of the 
self-energy vanishing within the gap region except for a 
narrow spike at the frequency where the real part of the 
self-energy changes sign. But the Z{0) parameter can 
still be defined on the imaginary axis! 

Conclusion We have discovered an approximate scal- 
ing relation for the electron-phonon problem that holds 
over a wide range of phonon frequencies, coupling 
strengths, and anharmonicities. It relates the wavefunc- 
tion renormalization parameter, extracted from the imag- 
inary axis, to the transition temperature. Our results 
show that anharmonic models that are tuned to the same 
Z{0) will show similar TcS as harmonic models with the 
same Z{0). We feel this underlying scaling behavior helps 
explain the success of the quasiharmonic approximation 
employed in describing nearly all real materials. 
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FIG. 1. Transition temperatures for the chessboard phase 
charge-density-wave order at half filling. (a) Tc for the 
harmonic case. The horizontal axis is the effective elec- 
tron-electron coupling \U\ — g^/ii- Three frequencies are in- 
cluded f7 = 0, SI = 0.5, and SI — oo. The dashed line is a 
conserving second-order perturbative approximation (includ- 
ing vertex corrections) and the solid line is a second-order 
strong-coupling approximation, (b) Tc for the anharmonic 
case. The horizontal axis is a measure of the anharmonic 
contribution to the energy (as described in the text). Open 
symbols are QMC simulations at f7 = 0.5, while solid symbols 
are exact solutions at f2 = 0. 



FIG. 2. Scaling curve for the anharmonic phonon prob- 
lem. The horizontal axis is the imaginary-axis-extrapolated 
wave-function renormalization parameter, extrapolated at Tc, 
while the vertical axis is the CDW Tc at half filling. In- 
cluded in the curve are weak-coupling results for the classical 
phonons with (dotted) vertex corrections, a conserving sec- 
ond-order approximation for the fl = oo case (dashed), and 
the approximate form described in the text (solid). 



FIG. 3. Interacting density of states for a generic 
weak-coupling harmonic case {Q = 0, g — 0.625 and 
Qan = 0.0) for six temperatures running from top to bot- 
tom at c^ = r = 0.04, 0.0625, 0.125, 0.25, 0.5, 1.0, and 2.0. 
The real and imaginary parts of the self energy are plotted in 
the insets. Notice how the slope of the real part of S has the 
wrong sign, and how the imaginary part has a local minimum 
at the chemical potential, signifying that this system is not a 
Fermi liquid. 



FIG. 4. Interacting density of states for a generic 
weak-coupling anharmonic case {il = 0, g — 2.5 and 
Oan ~ 4.5) for the same six temperatures shown in Figure 
3. Notice how the density of states is now asymmetric at high 
T but becomes more symmetric as T decreases. 



FIG. 5. Interacting density of states for a generic 
strong-coupling anharmonic case {fl = 0, g — 2.5 and 
(Xan = 0.004) for six temperatures running from top to bot- 
tom (at tj = 0) T = 0.075, 0.125, 0.25, 0.5, 1.0 and 2.0. 
Notice how the density of states develops a pseudogap at the 
chemical potential, which then becomes a true gap as T is 
lowered. 
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